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We study single electron transport across a single Bi dopant in a Silicon Nanotransistor to assess how the 
strong hyperfine coupling with the Bi nuclear spin 1 = 9/2 affects the transport characteristics of the device. 
In the sequential tunneling regime we find that at, temperatures in the range of lOOmif , dl/dV curves reflect 
the zero field hyperfine splitting as well as its evolution under an applied magnetic field. Our non-equilibrium 
quantum simulations show that nuclear spins can be partially polarized parallel or antiparallel to the electronic 
spin just tuning the applied bias. 



PACS: 73.23.Hk, 31.30.Gs, 74.55. +v, 75.75.-c 

The amazing progress both in the silicon process- 
ing technologies and in the miniaturization of silicon 
based transistors has reached the point where single- 
dopant transistors have been demonstrated. 1 7 Whereas 
this progress has been fueled by the development of clas- 
sical computing architectures, it might also be used for 
quantum computing. In this regard, the electronic and 
nuclear spins of single donors in silicon are very promis- 
ing building blocks for quantum computing. 8 10 Progress 
along this direction makes it necessary to implement sin- 
gle spin readout schemes both for electronic and nuclear 
spins. Single electronic spin readout has been demon- 
strated, both in GaAs quantum dots as well as in P doped 
Silicon Nanotransistors. 11,12 

The readout of the quantum state of a single nuclear 
spin, much more challenging, has been demonstrated for 
NV centers in diamond taking advantage of single spin 
optically detected magnetic resonance afforded by the ex- 
traordinary properties of that system. 13 Single nuclear 
spin readout with either optical 14 or a combined electro- 
optical techniques 15 has been proposed, but remains to 
be implemented. Here we explore the electrical readout 
of a single nuclear spin, more suitable for an indirect 
band-gap host like Si. A preliminary step is to construct 
a circuit whose transport is affected by the quantum state 
of the nuclear spin. There is ample experimental evidence 
of the mutual influence of many nuclear spins and trans- 
port electrons in III-V semiconductor quantum dots in 
the single electron transport regime. 16 19 In particular, 
Kobayashi et ai have reported hysteresis in the dl/dV 
upon application of magnetic fields, reflecting the real- 
ization of different ensemble of nuclear states coupled to 
the electronic spin via hyperfine coupling. 19 

Here we propose a device where a single nuclear spin 
is probed in single electron transport. We model the 
single electron transport in a silicon nanotransistor such 
that, in the active region, transport takes place through 
a single Bi dopant, see Fig. 1. We show that, at suffi- 
ciently low temperatures, the dl/dV curves of this de- 
vice probe the hyperfine structure of the dopant level. In 



FIG. 1. (color online) a) Scheme of the Si:Bi FinFET nan- 
otransistor. b) Trapping Coulomb potential of the Bi dopant 
an single energy level participating in the transport. 



turn, the occupations of the nuclear spin states are af- 
fected by the transport electrons. Whereas single dopant 
transistors have been demonstrated for single P, As and 
B, in Si, 3,4 ' 6,12 we choose Bi because it has a much 
larger hyperfine splitting, 20 22 due to both a larger nu- 
clear spin 7 = 9/2 and a larger hyperfine coupling con- 
stant (A « 6.1/ieV). The zero-field splitting of the Bi 
donor level is given by 5A and has been observed by 
electron spin resonance 20 22 and in photoluminescence 
experiments with many dopants. 23 

We consider the sequential transport regime, where the 
occupation of the donor level fluctuates between q = 
and q = 1. In the q = state, the nuclear spin interacts 
only with the external field. In the q = 1 state, the 
electron and the nuclear spin are hyperfine coupled. The 
Hamiltonian that describes both states reads 20 22,24 



e d + eV G + AS-I + huj e S z ) + hu N I z , (1) 



n 



where is the donor energy level with respect to the 
Fermi energy, which we take as Ep = 0, and Vq denotes 
an external gate voltage. We assume that valley degen- 
eracies of the donor level are split-off and neglect the 
valley degree of freedom. The third term is the hyper- 
fine coupling, and the last two, where huj e = g e ^BB z 
and hwN — 9nl^NB Zl correspond to the electron and 
nuclear Zeeman terms, with g e (g n ) the electron (nu- 
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FIG. 2. (color online) Scheme of the current-induced allowed 
transition for the a) q = 1 charged system and b) q — 
uncharged system. It has been assumed that Huon <C ksT <C 
^ e < A . 

clear) g- factors and /i# (/ijv) the Bohr (nuclear) mag- 
neton. In equilibrium, i.e., at zero bias, the occupation 
of the dopant level depends on the value of the addi- 
tion energy, which ignoring the Zeeman terms and the 
tiny correction due to the hyperfine coupling, is given by 
£o(V G ) = e d + eV G . 

We denote the q = eigenstates as \m). Their energies 
read as e m = hjj^I z . The eigenenergies and eigenvectors 
of q = 1 are denoted by cm and \M). The q = 1 zero- 
field Hamiltonian AI • 5 can be diagonalized in terms of 
the total angular operator F, resulting in two multiplets 
(F=4, F=5) with energies E F = 4 -11A/4 and E F = 5 
9A/4, and a zero-field splitting A = hA « 30/ieV. At 
finite magnetic field, the exact eigenvalues of 1-L can also 
be calculated analytically. 22 The corresponding energy 
levels are shown in Fig. 2. 

The tunneling Hamiltonian between the single Bi 
dopant level and the source and drain electrodes reads 
as 

Utun = ^V X {dlc X<T +h.c), (2) 
Act 

where operator c\ j(7 annihilates an electron with spin a 
and orbital quantum number A = 77, /c, with wave vec- 
tor k and electrode index 77 = S, D, while operator d a 
annihilates a spin a electron in the dopant level. The 
scattering rate for the tunneling process, ignoring the 
hyperfine coupling, is given by Tq = ^l^l 2 ^, where 
is the density of states of the electrode. Our model is 
very similar to the one used to describe single electron 
transport through a quantum dot exchanged coupled to 
a single Mn atom. 25 ' 26 

The dissipative dynamics of the electro-nuclear spin 
system under the influence of the coupling to the elec- 
trodes is described by a Bloch-Redfield (BR) master 
equation. 26 ' 27 The coupling to the reservoir, given by the 
tunneling Hamiltonian, involves transitions between the 



q = and q = 1 manifolds. The corresponding transition 
rates are be calculated using the Fermi golden rule with 
Htun as the perturbation: 26 

ri,M = rSEK M l J *( m WI 2 > (3) 

cr 

where \I Z , a) = \I Z ) (g) |cr). In the following we take the 
applied bias convention us — Pd = eV, with ps = eV/2 
and [Id = —eV/2. For a given temperature, bias and 
gate voltages and Hamiltonian parameters, we obtain the 
steady state solution of the master equation, ignoring the 
effect of the fast-decaying coherences. This yields the 
steady state occupations P m (V) and Pm(V). 

We consider the sequential tunneling regime, in which 
the energy level broadening induced by coupling to the 
electrodes To is small, hTo <C ksT . This also justifies the 
markovian approximation implicit in the Bloch-Redfield 
master equation. In this regime, current flows when the 
bias enables charge fluctuations of the dopant level. The 
steady state current corresponding to electrons flowing 
from the source electrode to the dopant level is given by 

I = eJ2 {P m (V)f s (A M , m )T s mM 

m,M 

-P M (V)[l-f s (A Mjrn )}T s mM }, (4) 

where A M ,m = e M - e m and f s (e) f(e - fi s ) is the 
Fermi function relative to the chemical potential of the S 
electrode. The first term in the right hand-side of Eq.(4) 
represents the electrons flowing from the S electrode to 
the empty Bi, while the second one corresponds to elec- 
trons flowing from the q = 1 Bi to the S electrode. In 
steady state, the continuity equation ensures that cur- 
rent between the dopant and the drain is the same than 
the source-dopant current. 

Figure 3a) shows the differential conductance dJ[/dV 
map for zero- applied magnetic field, with / = I/(^o) 
and Tq = To/2. At zero bias, the conductance is zero 
except at the special value of Vq for which the addition 
energy vanishes. Far from this point, the zero-bias charge 
of the dopant state, hereafter denoted with go, is either 
qo = or qo = 1. The finite bias conductance has a 
peak whenever the bias energy, eV/2, matches the energy 
difference between two states with different charge, m for 
q = and M for q = 1, that are permitted by the spin 
selection rule implicit in Eq.(3). The height of the peak is 
proportional to both the non-equilibrium occupations P m 
and Pm and to the quantum mechanical matrix element 
M' This determines the very different spectra when 
the zero bias charge in the dopant is q = or q = 1 . The 
width of the dl/dV peaks is proportional to feT, so that 
the dl/dV spectra can resolve the hyperfine structure 
provided that ksT is smaller than the splitting of the 
levels. The energy differences inside the F = 4 and F = 5 
manifolds, see Fig. 2), are roughly proportional to A. 
Thus, while the zero-field splitting can be resolved at 
T = 0.3 K, temperature must be significantly below 50 
mK to resolve the finite field structure, see Fig. 3c). 
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FIG. 3. (color online) a) and b) Contour plot of the d^/dV 
vs. applied bias V and on-site energy so at zero magnetic field 
(left) and B — 0.6T (right), c) and d) Conduction spectrum 
dI_/dV as a function of applied bias for different magnetic 
fields at so = — 0.4/xeV and so = 0.8meV respectively. White 
horizontal lines in panel a) and b) marks the values of So in the 
2D plots c) and d). In all cases, T = lOmK and hVo = 0.1/xeV. 



Let us consider first the qo = 1 case (left panel in 
Fig. 3). At 10 mK only the ground state(s) is (are) oc- 
cupied. Thus, a single transition is seen, from the q = 1 
to the q = states. As the magnetic field is ramped, 
the energy of the transition increases, reflecting the elec- 
tronic Zeeman shift. In contrast, in the qo = case (right 
panel in Fig. 3), all the Zeeman split nuclear levels are 
equally populated, even down to mK temperatures. Spin 
conservation selection rule implicit in Eq. (2) connects 
these 10 quasi-degenerate states of the q = manifold to 
the hyperfine spin-split levels of the q = 1 manifold with 
different energies. As a result, the dl/dV curve reveals 
2 peaks at zero field, reflecting the splitting between the 
F = 4 and F = 5 states. At higher fields, the two zero- 
field peaks split in up to 10 peaks, that can be resolved 
at low enough temperature [see Figs. 3c) and 4b)]. 

Interestingly, the application of a bias to the qo = 
state, for which the nuclear spin states are randomized, 
can result in a finite average nuclear magnetic moment. 
We show this in Fig. 4a) for finite B. At zero bias, the 
charge of the dopant level is qo = 0, and the nuclear spins 
are randomized. When the bias hits the addition energy 
a selective depopulation of a given I z level of the q = 
manifold starts, in favor of a q = 1 state that mixes the 
I z and 4±1 components, resulting in a net accumulation 
of nuclear spin. When all the transitions to the F = 4 
manifold are allowed, the nuclear spin vanishes again. 
Then, when the bias permits the transitions to the F = 
5 manifold, the nuclear spin accumulation starts in the 
opposite direction. Thus, when \eV/2\ matches the center 
of the F = 4 multiplet, see Fig. 4a), the nuclear spins 
tend to align antiparallel to the electronic spin. Then, 
when \eV/2\ reaches the center of the F = 5 multiplet, 
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FIG. 4. (color online) a) Average electronic occupation of 
the Bi, (Q)/\e\ (black line) and nuclear and electronic spins, 
(I z ) (red line) and (S z ) (blue line), respectively, b) d^/dV 
vs. bias for different temperatures. Same parameters as 3c) 
with B 0.6T. 



the nuclear spins prefer aligning parallel to the electronic 
spin. 

Whereas all our results discussed so far refer to steady 
state conditions, it is worth pointing out that there are 
two very different time scales in the dynamics of the sys- 
tem. Whereas the charge equilibrates in the dopant level 
in a time scale set by I/Tq, the nuclear spin relaxation, 
dominated by many events of hyperfine exchange with 
the electronic spin and subsequent recharging of the Bi, 28 
takes place at a much longer time scale, hundreds of time 
larger than 1/To, but still much shorter than the intrinsic 
T\ time of the nuclear spin. Thus, charge fluctuations in 
the Bi induce nuclear spin relaxation. 28 

We finally discuss the experimental feasibility of our 
proposal with state of the art techniques. First, accord- 
ing to our simulations, see Fig. 4b), the finite field hyper- 
fine splitting is resolved at 10 mK but not a 20 mK. At 40 
mK the 2 humps associated to the F = 4 and F = 5 man- 
ifolds are clearly resolved. Keeping the transport in the 
sequential tunneling regime requires that hTo <C fc^T, 
which at lOmK, translates into / <C 200pA. This is within 
reach of experimental setups. 12,16,19,29,30 

In conclusion, we have studied the single electron 
transport spectroscopy of the hyperfine structure of a Bi 
dopant in a silicon nanotransistor. We have shown that, 
at sufficiently low temperatures, and when the dopant 
is ionized with a gate, the dl/dV corresponding to se- 
quential transport can resolve the hyperfine spectrum of 
the electron in the donor level. In addition, the non- 
equilibrium transport at finite field results in a hyper 
polarization of the nuclear spin state, or nuclear spin ac- 
cumulation. These results are different from our previous 
work, where we considered the same system in a different 
transport regime, cotunneling, and we showed that in- 
elastic cotunneling of the dopant in the q = 1 state could 



also resolve the hyperfine spectrum and drive the nuclear 
spin states out of equilibrium. 24 Future work should de- 
termine how, in the cotunneling regime, the appearance 
of the Kondo effect 7 competes with the reported effect. 
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